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Abstract. The spiral is one of Nature’s more ubiquitous shape: it can be seen in 
various media, from galactic geometry to cardiac tissue. In the literature, very specihc 
models are used to explain some of the observed incarnations of these dynamic entities. 
Barkley [1,2] first noticed that the range of possible spiral behaviour is caused by the 
Euclidean symmetry that these models possess. 

In experiments however, the physical domain is never perfectly Euclidean. The 
heart, for instance, is hnite, anisotropic and littered with inhomogeneities. To capture 
this loss of symmetry (and as a result model the physical situation with a higher degree 
of accuracy), LeBlanc and Wulff introduced forced Euclidean symmetry-breaking 
(FESB) in the analysis, via two basic types of perturbations: translational symmetry¬ 
breaking (TSB) and rotational symmetry-breaking terms. In [3,4], they show that 
phenomena such as anchoring and quasi-periodic meandering can be explained by 
combining Barkley’s insight with FESB. 

In this article, we provide a fuller characterization of spiral anchoring by studying 
the effects of n simultaneous TSB perturbations, where n > 1. 
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1. Introduction 

Spiral waves have been observed in a variety of experimental contexts, ranging from 
the well-known Belonsov-Zhabotinsky chemical reaction to the electrical potential in 
cardiac tissne [1,2,5-16]. In this last case, spiral waves are believed to be a precnrsor to 
several fatal cardiac arrythmias {e.g. ventricular tachycardia and ventricular hbrillation) 
[13,17,18]. A thorough understanding of the various dynamical properties of spiral waves 
is therefore warranted. 

One of the most interesting and fruitful approaches in recent years to the study 
of spiral waves has been to use the theory of equivariant dynamical systems to derive 
hnite-dimensional models for many of the observed dynamical states and bifurcations of 
spirals. The pioneer of this approach was Barkley, who realized that the experimentally- 
observed transition from rigid rotation to quasi-periodic meandering and drifting 
could be explained using only the underlying symmetries (the group §E(2) of all 
planar translations and rotations) of the governing reaction-diffusion partial differential 
equations: he derived an ad hoc system of 5 ordinary differential equations with SE(2) 
symmetry which model a Hopf bifurcation from a rotating wave, and then showed 
that this hnite-dimensional system replicated the experimentally-observed transition to 
meandering and drifting [1,2,19]. Sandstede, Scheel and Wulff later proved a general 
center manifold reduction theorem for relative equilibria and relative periodic solutions 
in spatially extended inhnite-dimensional SE(2)-equivariant dynamical systems, thereby 
providing mathematical justihcation for Barkley’s approach [20-24]. 

One of the advantages of this equivariant dynamical systems approach is that 
one can often give universal, model-independent explanations of many of the observed 
dynamics and bifurcations of spiral waves. For example, the above-mentioned Hopf 
bifurcation from rigid rotation to quasi-periodic meandering and drifting has been 
observed in both numerical simulations [16] and in actual chemical reactions [5]. Another 
example is the anchoring/repelling of spiral waves on/from a site of inhomogeneity, 
which has been observed in numerical integrations of an Oregonator system [8], in 
photo-sensitive chemical reactions [7] and in cardiac tissue [11]. Using a model- 
independent approach based on forced symmetry-breaking, LeBlanc and Wulff showed 
that anchoring/repelling of rotating waves is a generic property of systems in which 
the translation symmetry of SE(2) is broken by a small perturbation [3]. Similarly, 
some dynamics of spiral waves observed in anisotropic media {e.g. phase-locking and/or 
linear drifting of meandering spiral waves) have been shown to be generic consequences 
of rotational symmetry-breaking [4,12,25-28]. 

Consider as a paradigm a system of reaction-diffusion partial differential equations 


du 

'dt 


D-V‘^u + f{u) 


( 1 . 1 ) 
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where m is a fc-vector valued function of time and two-dimensional space, D is a 
matrix of diffusion coefficients and / : is a smooth reaction term. Many 

of the phenomena in which spiral waves are observed experimentally are modeled by 
systems of the form (HD). Moreover, Scheel has proved that systems of this form can 
admit time-periodic, rigidly rotating spiral wave solutions [29]. Implicit in the form of 
equations (HH) is the fact that the medium of propagation is completely homogeneous 
and isotropic. Mathematically, this is represented by the invariance of (EH) under the 
transformations 

u(t, x) I—>• u(t, xi cos 9 — X2 sin 9 + pi,Xi sin 9 + X2 cos 9 -|- P2), ( 1 . 2 ) 

where {9,pi,p2) G x and x G [30,31]. The group of all transformations of the 
form m is isomorphic to the special Euclidean group SE(2) of all planar translations 
and rotations. 

When studying the effects of inhomogeneities on the propagation of spiral waves, 
one must consider a larger class of models than (EH), since inhomogeneous media do 
not possess Euclidean invariance. For example, one might consider systems of the form 

— = D - W^u + f{u) -F Xg{u, ||x|p. A) (1.3) 

which are perturbations of EH)- Such systems could model a spatially extended 
react ion-diffusion medium in which there is one site of inhomogeneity (with circular 
symmetry) centered at the origin of For instance, the Oregonator model which 
is used to study spiral anchoring in [ 8 ] is of the form (ll.3|l . When A 7 ^ 0, (II. 3|) has 
rotational symmetry about the origin, but does not possess any translation symmetry. 
This phenomenon is called forced translational symmetry-breaking] it is studied in detail 
in [3]. 

In this paper, we use a similar equivariant dynamical systems approach to study 
the problem of spiral wave dynamics (specihcally, with regards to anchoring/repelling) 
in media in which there are several sites of inhomogeneities (as opposed to just one site), 
of which cardiac tissue is an important example. 

We will make several simplifying assumptions which are meant to make the analysis 
more tractable. First, we assume that the inhomogeneities consist of a hnite number of 
“sources” which are localized near distinct sites Ci,..., ^^ in the plane. Second, we will 
assume that these n sources of inhomogeneity are independent in the following sense: 
we introduce n independent real parameters Ai,..., A„ which give some measure of the 
relative “amplitudes” of the sources. In particular, when all the Xi are zero except, say 
Aj* 7 ^ 0 , then there is only one source of inhomogeneity localized near the point 
In that case, we will make a third simplifying assumption: the single inhomogeneity 
is circularly symmetric around the point Ci*- The following is an example of a class 
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of reaction-diffusion partial differential equations which are perturbations of (mi) and 
which might model such a situation; 


ut = DAu + f{u) + ^ \j 


Dji\\x - Cjf, \)Au + fj{u, ||x - Cjf, A) 


(1.4) 


where the functions Dj, fj are bounded and smooth enough. The goal of this paper is to 
provide a detailed analysis of a larger class of abstract dynamical systems which share 
the symmetry properties of dEl: 

(51) when Ai = ■ ■ ■ = A^ = 0, the systems are invariant under the action (j1 .‘2fl of the 
group SE(2), 

(52) when all the Aj are zero except Aj*, the systems have rotational symmetry about 
the point but they do not generically possess translation symmetries, 

(53) when two or more of the Aj are non-zero, the systems do not generically possess 
any of the symmetries dm except for the identity. 

Our results will apply to the subclass LCq of systems whose members also generate a 
smooth local semi-flow on a suitable function space [20-24], as well as some technical 
conditions which will be specihed as we proceed. 

The paper is organized as follows. In the second section, we derive the center bundle 
equations of the semi-flow of a system in LCq, near a hyperbolic rotating wave. We state 
and prove our main results in the third section: to wit, spiral anchoring is generic in 
a parameter wedge. Then, we provide a visual criterion characterizing the anchoring 
wedges in the case n = 2. Finally, we perform numerical experiments demonstrating 
the validity of our results. 


2. Reduction to the Center Bnndle Equations 

Let X be a Banach space, U C M” a neighborhood of the origin and be a smoothly 
parameterized family (parameterized by A G W) of smooth local semi-flows on X. 

Let SE(2) = C-i-§0(2) denote the group of all planar translations and rotations, 
and let 

a : §E(2) —^ GL(X) (2.1) 

be a faithful and isometric representation of SE(2) in the space of bounded, invertible 
linear operators on X. For example, if AC is a space of functions with planar domain, a 
typical SE(2) action (such as ()1.2j) in the preceding section) is given by 


(a(7)M)(x) = m(7 ^(x)), 7 G §E(2). 
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We will parameterize SE(2) as follows: SE(2) = C x S^, with multiplication given by 
{Pi, Pi) ■ {P 2 , P 2 ) = + Pi, Pi + <^ 2 ), V (pi, ipi), {p 2 , P 2 ) ^ C X Sh For fixed ^ G C, 

we define the following subgroup of SE(2): 

§o( 2 )^ = {(e,o)-(o, 0 )-(-e,o) 1 

which is isomorphic to SO(2), and represents rotations about the point We will 
assume the following symmetry conditions on the family of semi-flows. 

Hypothesis 1 There exists n distinct points ^i,... ,^n in C such that if ej denotes the 
vector of the canonical basis in M", then Vn G X, a 7 ^ 0 , f > 0 , 

7 G SO(2)^., and 
Vy G §E(2). 

Hypothesis [D basically states that (a) when A = 0, the semi-flow <Fi,o is SE(2)- 
equivariant; (b) when A 7 ^ 0 is near the origin and along the coordinate axis of M"’, 
the semi-flow is only SO(2)^^-equivariant (i.e. it only commutes with rotations about 
the point ^j), and (c) when A is not as in (a) or (b), the semi-flow has (generically) 
trivial equivariance. 

We are interested in the effects of the forced symmetry-breaking on normally 
hyperbolic rotating waves. Therefore, we will assume the following hypothesis. 

Hypothesis 2 There exists n* G X and Tl* in the Lie algebra of §E(2) such that 
is a rotation and = a{e^*^)u* for all t. We also assume that the set 

{ A G C I |A| > 1} zs a spectral set for the linearization a{e~^*)D^i Q{u*) with projection 
F* such that the generalized eigenspace range(F,,) is three dimensional. 

For sake of simplicity, we will only be interested in one-armed spiral waves; there¬ 
fore, we assume that u* in hypothesis El has trivial isotropy subgroup. While hypotheses 
1 and 2 hold for a large variety of spirals (such as decaying spirals), there is also a large 
family of spirals for which they don’t (including Archimedean spirals) [29].§ 

Let LCo be the collection of all abstract dynamical systems that do satisfy them, as well 
as all other hypotheses required in order for the center manifold theorems of [20-23] to 
hold, and let be produced by some member of LCq. It follows that for A near the 
origin in M”, the essential dynamics of the semi-flow near the rotating wave reduces 

§ However, even in the case of Archimedean spirals (for which hypothesis 1 fails), finite-dimensional 
center-bundle equations which share the symmetries of the underlying abstract dynamical systems have 
been shown to possess a definite predictive value in terms of possible dynamics and bifurcations of these 
spiral waves [1-5]. 
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to the following ordinary differential equations on the bundle C x (see [32] for more 
details): 

P = e^^V + GP{p,p,Lp,\) 
p = uj + G'^{p,p,ip,X) 

where Id is a complex constant, ci; 7 ^ 0 is a real constant, G^ and G‘^ are smooth, 
uniformly bounded in p, and such that GP{p,p, <p, 0) = 0 and G'^{p,p, <p, 0) = 0. If A is 
near the origin, we can re-scale time along orbits of (EH to get 

p = e*^w + ^(p,p,(p,A) 
p = 1 

where Q is smooth, uniformly bounded in p, and such that G{p,p, <p, 0) = 0. Of course, 
Q is not completely arbitrary because of the symmetry conditions in hypothesis ^ A 
simple computation and Taylor’s theorem lead to the following. 


Proposition 2.1 The symmetry conditions in hypothesis^ imply that the equations 
(Q have the general form 


p 




(p - A) 

i=i 


(2.4) 


where, without loss of generality, (p(t) = t, v E C, X = {\i,, Xn), and the functions 
H, are smooth and uniformly bounded in p. 

A 27r—periodic solution px of (Q is called a perturbed rotating wave of (j2.4j) . Dehne 
the average value 

1 

[Pa]a = ^ (2.5) 

If the Floquet multipliers of p\ all he within (resp. outside) the unit circle, we shall say 
that [pa]a is the anchoring (resp. repelling, or unstable anchoring) center of px. 

In the following section, we will perform an analysis of anchoring of perturbed 
rotating waves of (Q for parameter values near A = 0. 


3. Analysis of the Center Bundle Equations 

Equations (Q represent the dynamics near a normally hyperbolic rotating wave 
for a parameterized family ^t x of semi-flows satisfying the forced-symmetry breaking 
conditions in hypothesis Cl We start with a brief review of the case n = 1 which was 
studied in detail in [3], and then present new results on the general n case. 
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In this case, we may assume without loss of generality that = 0, so that dZl has the 
form 


P = e'* [v + XH{pe-^\pC\ A)] (3.1) 

where A G M is small. By writing w = iv, this system becomes 

w = —iw + XH{w,w, X) (3.2) 

where H{w,w, A) = H{w — iv,w + iv, A). The following theorem is proved in [3]. 

Theorem 3.1 Let a = Re(iAii/(0, 0, 0)), where H is as in If a ^ 0, then for all 

A 7 ^ 0 small enough, <\3. ijl has a hyperbolic rotating wave 

p(f) = (-in + 0(A)) e**, (p(t)=t. (3.3) 

The origin [p]a = 0 is an anchoring center if aX < 0; it is a repelling center if aX > 0. 

Remark 3.2 In the case where the semi-flow ^ is generated by a system of planar 
react ion-diffusion partial differential equations, the solution (Q represents a wave 
which is rigidly and uniformly rotating around the origin in the plane. In the case where 
aA < 0, the rotating wave is locally asymptotically stable. When aA > 0, the rotating 
wave is unstable (see [8] for an experimental characterization of this phenomenon in an 
Oregonator model). 

3.2. The Case n > 1 

One might think that the combination of many perturbations would just combine the 
effects of each perturbation, so that spirals would be observed anchoring at each of the 
centers, but we shall see that this is not usually the case. 

By re-labeling the indices in (jZl if necessary, we can temporarily shift our point 
of view so that plays the central role in the following analysis. Then, under the 
co-rotating frame of reference z = p — fi + ie^^v, dZl becomes 

n 

z = p — = e** ^jHj [{z — iv, {z — Ci)e**+ iT, A), (3.4) 

where (j = for j = 1,..., n. 

When Ai 7 ^ 0 and A 2 = ■ ■ ■ = A^ = 0, we find ourselves in the situation described 
in the previous subsection. Now, set e = Ai, = 1 and Xj = pjS for j = 2,... ,n and 



Spiral anchoring in media with multiple inhomogeneities 


p = (yU. 2 ,..., Pn) G M” Then (jH.4|l can be viewed as a perturbation of the correspond¬ 
ing equation in the case n = 1. Note that Ci = 0 and A = ( 1 , p)e. 


Equation (Q rewrites as 


n 

z = ee** ^ pjHj {{z - 0)e“* - iv, {z - 0)e*V iv, (1, p)e). (3.5) 

Let Hj{w, w, e, p) = Hj (w — iv,w + iv, (1, p)e) for j = 1,.. .n. Then ()3.5|1 becomes 

z = ee^^K {ze~^^, ze**, t, e, p) (3.6) 


where K{w, w, t, e, p) = pjHj{w — Qe w — s, p) is 27r—periodic in t. 

i=i 

Set = DiHi{—iv,iv,Gi). The time—27r map P of (I3.6|l is given by 


P{z, z,e,p) = z + 27re oiz + 0(\z\^) + 0[e, p 2 ,..., Pn) 


(3.7) 


near z = 0 and (e, p) = ( 0 , 0 ). 

Hyperbolic hxed points of (EH) correspond to hyperbolic 27r—periodic solutions 
of (lT6ll . and so to perturbed rotating waves of (El, that is, the path traced by the 
solution wave need not be circular. As z = 0 is not generally a hxed point of dH, these 
perturbed rotating waves may not be centered at ^i. Indeed, let 


B{z, z, £, p) = aiz 0{\z\^) -1- 0(£, p 2 , • • •, l^n) (3.8) 

be the function inside the square brackets in (13.7|) . Note that B{0, 0, 0, 0) = 0 and that, 
generically, DiB{0, 0, 0, 0) = ai 7 ^ 0. By the implicit function theorem, there is a unique 
smooth function z{e, p) dehned near {e, p) = ( 0 , 0 ) with z{0, 0 ) = 0 and 

B{z{e,p),z{e,p),e,p) =0 (3.9) 

near ^ = 0. This leads to the following theorem. 

Theorem 3.3 Let ai he as in the preceding discussion, with Re(ai) 7 ^ 0. If the 
parameters are small enough to satisfy the conditions outlined in the proof below, the 
time—2% map has a unique family of hyperbolic fixed points, whose stability is 

exactly determined by the sign ofeRe{ai). 

Proof: Let B be as in dH and z{e,p) be the unique continuous function solving the 
equation H = 0 for small parameter values, as asserted above. When £ = 0, any point 
of is a non-hyperbolic hxed point of P and so, from now on, we will assume that 
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e 7 ^ 0. If that is the case, and if £ and ||/i|| are small enough, the eigenvalues 
of DP{z{e, p),e, fi) satisfy 

= 1 + 47r£Re(ai) +eO{e,fi) ^ 1, 

since Re(ai) ^ 0. In other words, the fixed point z{e,fi) is hyperbolic. When 
eRe(ai) < 0 , the eigenvalues lie inside the unit circle and the hxed point is asymp¬ 
totically stable; otherwise, it is unstable. □ 

We are now able to formulate and prove the following result. 

Theorem 3.4 Suppose the hypotheses of theorem AS.,ti are satisfied. Then there exists a 
wedge-shaped region near X = 0 of the form 

Wi = {(Ai,..., Xn) G M"' : |Aj| < lRij|Ai|, Wij > 0, for j ^ 1 and Ai near 0 } 

such that for all 0 ^ X E Wi, has a unique perturbed rotating wave with 

center [iSj[]A generically away from fi. Furthermore, [iS{]a is a center of anchoring 
when Ai Re(ai) < 0. 

Proof: For j 7 ^ 1, let hFij > 0 be such that the conclusion of theorem 13.31 holds for 
any pj with \jjij\ < IFij. Let Wi be as stated in the hypothesis. If (e,//) is such that 
the time—27r map (j3.7j) has a hyperbolic fixed point z{e,p) near 0, then (I3.5|l has a 
hyperbolic 27r—periodic orbit centered at a point near z = (i = 0. 

For j 7 ^ 1, let Ai = e 7 ^ 0 be small enough and set Xj = fijS. Then A G Wi, as 

\Xj\ = \iaj\\e\ < IFijIAil for j 7 ^ 1, 

and z^^fit) is a 27r—periodic orbit for the parameter A, which we denote by zx(t). Since 
p = z — ie^^v + fi, (El has a unique perturbed rotating wave with 

[^IIa = ^ j + ^1) dt = ^i + [^a]a. 

If 0 7 ^ A G W is such that pj = Aj/e 7 ^ 0 is fixed for j = 2,..., n, then [ 2 :a]a = 0(1) as 
Ai —*■ 0 and so [iS1]a 7 ^ ■Ci) generically. The conclusion about the stability of follows 
directly from theorem 13.31 □ 

Remark 3.5 When A approaches the Ai—axis away from the origin, [iS{]a —^ On 

the other hand, when the parameter values stray outside of Wi, all that can generically 
be said with certainty is that solutions of dH locally drift away from ,^ 1 , which cannot 
then be a center of anchoring. After drifting, the spiral may very well get anchored at 
some point far from ,^ 1 , depending on the global nature of the perturbation functions 
Hj in ()2.4|) . 
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The preceding results have been achieved by considering (El under a co-rot at ing frame 
of reference around Of course, since the choice for ,^i was arbitrary, corresponding 
results must also be achieved, in exactly the same manner, when the viewpoint shifts 
to another For j = 1,... ,n, let aj = DiHj{—iv, iv, 0) be the anchoring coefficients 

of El . 

Theorem 3.6 Let k G // Re(afc) ^ 0, then there exists a wedge-shaped 

region near X = 0 of the form 

Wk = {(Ai,..., An) G M" : \Xj\ < Wk,j\Xk\, Wk,j > 0, for j ^ k and Xk near 0 } 

such that for all 0 X E Wk, (H has a unique perturbed rotating wave with 

center [5^]a generically away from fk- Furthermore, [5^]a Is a center of anchoring 
when AfcRe(afc) < 0. 

Clearly, the remark that appears after the proof of theorem 1.3,41 still holds. 

4. Characterization of Spiral Anchoring (n = 2) 

In the previous section, we described the (local) behaviour of spiral anchoring in small 
wedges around the parameter coordinate axes. In this section, we present a fuller char¬ 
acterization of spiral wave anchoring for the case n = 2.|| 

Let 0 ^ G M^, Aq = (Ai,0),A^ = (0, A 2 ) G and let P : x be 

a real analytic map with P{x, 0) = x, DP{x, 0) = I 2 for all x G R^, satisfying the 

following conditions: for rj G {0,^}, 

(PI) 3a;* > 0 such that P{ri,Ari) = 0, for all ||A^|| < a;*; 

(P2) the eigenvalues of DPffi, A^) lie both outside or both inside the unit circle for all 

0 ^ I I-A-t?! I ^ 

(P3) there is a wedge region surrounding the coordinate axis generated by A^ in 
parameter space (see figure 1) in which P has a (locally) unique manifold a:^(A) 
such that, for all A G w,^, 

(a) P(a;^(A), A) = x^(A); 

(b) Xn{X) —> ?7 as A approaches the coordinate axis away from the origin; 

(c) x^(A) shares its stability with p in (P2). 

When the hypotheses of theorem l3.6l hold. the associated time—27r map ()3.7|1 (viewed in 
real coordinates) satisfies (PI) —(P3). Numerous questions cannot be answered by local 
analysis alone. For instance: 


II Most of the analysis can be extended and adapted to the general case n > 2, but at the cost of 
substantial algebraic complications. 
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Figure 1. Wedges in parameter space corresponding to property (P3). 


(i) Can the wedges overlap? What does that imply for anchoring in (E3D? 

(ii) Can a wedge contain its “opposite” coordinate axis? 

(iii) If the wedges do not overlap, what is the nature of their complement? 

(iv) If there is a complement with non-trivial measure, what kinds of dynamics can be 
expected as the parameter vector A traces a circle around the origin in parameter 
space? 

We will provide answers to these questions by hrst studying a specihc map, then 
extending our results to the general mapping. 

4.I. A Specific Mapping 

Consider the mapping P : x given by 

P{x, X) = x + 27r[AiPo(a;) + ^ 2 G^{x)] , (4.1) 

where 0 7 ^ G M^, and Fq, are real analytic functions of x, A G 

Such a map is obtained by truncating the A—terms of order > 2 from the time—27r 
map ()d.7j) . for instance. According to theorem Id.hi the jacobians DFq{0) and DG^{^) 
have a particular structure. 

Proposition 4.1 If Fo{0) = 0, = 0? */ 

BF„(0) = -fj ani DGdO = 

where a,c ^ 0, then there exists a;* > 0 such that the map defined by o satisfies the 
conditions (PI) —(P3). 
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A{x) 


Fo(x'] Gj(x) 


(4.2) 


Then, x is a fixed point of iH) for A G if and only if A(x) ■ A = 0 , that is if and 
only if A G = ker 74 (x). Let (x, A) be such a pair. According to the implicit function 
theorem, as long as 

det A) - /) = dTT^ det (^AiBFo(x) + A 2 BG^(x)^ ^ 0, (4.3) 


there is a neighbourhood IT of A and a unique analytic function X : IT —> such that 
X(A) = X and A(X(A)) ■ A = 0 for all A G IT. By construction, X(A) is a hxed point of 
(j4.1|l for all A G IT. 

If diniL^ = 0 as a manifold, then = {0}. Consequently, the preceding implicit 
function theorem construction fails, which contradicts property (P3). We need thus only 
investigate hxed points x for which diniL^ ^ 0 . As the quantities under consideration 
are analytic, it can further be assumed that rankA(x) = 1 and diniL^ = 1 . 


We now show how to optimally extend the wedge regions using property (P3). 
Let (x*. A*), (x*. A*) G X (M^ — {0}) be such that x*, x* are hxed points of ()4.1|1 . 
A* G La;*, A* G La;*, and dOl) is satished for both pairs. According to the implicit 
function theorem, there are open neighbourhoods IT*, IT* of A*, A* G respectively, 
and a pair of unique real analytic functions X* : IT* —>■ X* : IT* —>■ for which 

X*(A*) = X*, X*(A*) = X* and 

A(X*(A)) ■ A = 0, for all A G W*, A(X*(A)) ■ A = 0, for all A G W*. 


Lemma 4.2 If A* G IT* = IT* fl IT* is such that X*(A*) = X*(A*), then X* = X* 

on it; 

Proof: The assertion follows from the uniqueness of the real analytic functions X*, X* 
in the implicit function theorem. □ 

Denote the punctured open disc of radius a;* centered at the origin by 5(0,u;*). Let 
?7 G {0, cn* > 0 be as in proposition 14.II and 0 7 ^ A,^ G be a point on the appropri¬ 
ate coordinate axis, as in properties (PI) and (P2). According to these same properties, 
?7 is a hxed point of (ED for A^j and det {D^P^rj, A^) — J) 7 ^ 0 . 

Lemma then implies the existence of a maximal open region W,^, dehned as a 
union of open sets IT C 5(0, cn*) (in much the same way as the maximal interval is 
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Figure 2. Wedge angles, with optimal wedge-like regions in parameter space. 


built in the Fundamental Theorem of ODE [33]), containing 0 5(0, a;*) and for which 
there is a unique real analytic function X,, : satisfying Xn = X^|w,,, where Xr^ 

is as in property (P3). Since x is a fixed point of (I4.1|l for 0 7 ^ A whenever A G L^, \Nn 
is described (in polar coordinates) 

W, = {(r,0):O< 

W, = {(r,0):O< 

where G ( 0 , 7 r/ 2 ] and 


In the latter case, we will say that 
titles G ( 0 , 7 r/ 2 ] are called the fore-angle and post-angle of IF^, respectively (see 

figure 2 ). 

The implicit function theorem fails to extend A{Xrj{X)) ■ A = 0 (that is, it fails to 
extend W,,) at {x*, X*) if either 

(Cl) det {XlDFo{x*) + X^DG^^x*)) = 0 and X,,(A) —> x* as A —*■ A*, or 
(C 2 ) ||X^(A)|| ^ 00 as A ^ AO 

Such events will be referred to as fold and 00 —catastrophes, respectively, or catastrophes, 
collectively. 

Let 0 < p < cu* and set 


by either one of 

r < a;* and < 6* < -h } 

r < O'* and 9 G [0, 27r]} 


0 if ?7 = 0 , 
jr /2 = 


(4,4) 


is catastrophe-free. In the former case, the quan- 


44(d =p(cos(s),sm(s))^. 


(4.5) 
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Assume is not catastrophe-free. Starting at (p, s^) G W^, denote the angles in 
( 0 , 7 r/ 2 ] measuring the first clockwise and the first counter-clockwise occurrence of a 
catastrophe along ■jp by 9~ and 9^ respectively. Then, ± 9^. 


4 . 1 . 2 . Fold bifurcation points. Modulo a simple regularity condition (see below), (Cl) 
is equivalent to the existence of a fold bifurcation curve in parameter space for iH). 
Indeed, in that case, (x*. A*) is a solution of 

P{x, X) — X = 0, det {DxP{x, A) — /) = 0. (4.6) 

If the (full) Jacobian of the left-hand side of (j4.6jl has rank 3 at that point, (USD has a 
fold bifurcation curve through A* [34]. Such solutions are in one-to-one correspondence 
with regular solutions of 


A{x) ■ A = 0, det {Dx [A(x) ■ A]) = 0. 


(4.7) 


Set 


ho — 


1 0 
0 0 


hi — 


0 0 
0 1 


I = 


0 1 
1 0 


ei = 




and 62 = 


0 


Cl, H 2 {x)= hoh{x) + hih{x)I 


and define TTi, ^ by 

Hi{x) = /loA(x) -h /oiA(x)J 
A quick computation shows that (HID can be written as 

A{x) ■ A = 0 , A''^Q(x)A = 0 , 


62- 


(4.8) 


where 


and 


Q{x) 


B{x) |C'(x)A 
lC{x) E{x) j 


B{x) = det DFo{x) 

C{x) = det DHi{x) + det DH 2 {x) 
E{x) = det DGhh- 


(4.9) 


Let X* be a fixed point of (SH) and denote = {A : A^(5(x)A = 0}. Generically, K^* 
consists of a single line or a pair of intersecting lines through the origin in parameter 
space. Writing Lj. = fl 5(0,cn*) and K^, = H 5(0,u;*), we can summarize the 
situation with the following proposition. 


Proposition 4.3 If {x*,X*) is a regular solution of o with {0} h Lx* ^ Kj,*, then 
(x*. A) is a fold bifurcation point of o for all A G L 3 ;*. 
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4-2. The Visual Criterion 

Set 3 = {( 2 ^, A) : P{x, X) = X and A G L 3 ; 7 ^ {0}} and 

fi:(3) = {x : 3X ^ 0 such that {x, A) G 3}- 

By construction, fi;(3) is the zero-set of dety4(x) in and 0,^ G fi:(3)- Generically, 
«:(3) is a collection C of isolated planar curves, whose constituents come in two vari¬ 
eties: bounded or unbounded.^ Denote this partition by C = Cb U Coo and let Co, C^ be 
the curves in C for which 0 G Gq and ^ & C^. 

Let 7 p : [0, 27r] —be the circle of radius p around the origin, parameterized as 
in iH). For each (x. A) G 3, define and as the intersection of that circle with 
and Kj;, respectively, and let Pp -.M? x [ 0 , 27r] —»■ be given by 

Pp{x,s) = a;-|-27rp[cos(s)Fo(a;) 4 -sin(s)Gj(a;)]. (4-10) 

Then consists of two antipodal points {Pa^^p}, and the hxed points (x, s) of Pp are in 
one-to-one correspondence with the ‘lines’ of hxed points [x, L^,) of P for which Lj, 7 ^ {0} 
(see proposition WM- 

Set 3p = {(a^, s) G X [0, 27r] : Pp{x, s) = x and 7 p(s) G Lg} and 

Kp(3p) = {a: : (x, s) G 3p}- 

By construction, Kp(3p) = «^(3)- Thus for each G G C, k~^{C) is a branch of hxed points 
in the bifurcation diagram of Pp. According to section 14.1.11 the converse also holds: 
each branch of hxed points in the bifurcation diagram of Pp projects down via Kp to a 
curve in C. 

The existence and location of fold catastrophes cannot be read directly from C, but 
the next proposition remedies that situation. 

Let {x*,a) G X (M^ — {0}) be such that ||q;|| = p, a G L^* C K^*. Recall that 
A{x*) 7 ^ 0. Then, {Aj^i{x*) Aj^ 2 {x*)) 7 ^ 0 for some j G {1,2} and 

Lx* = {A : Aj^i(x*)Xi -|- Aj^ 2 {x*)X 2 = 0}. 

The function Tj : dehned by 

rj(a;) = [Al2{x)B{x) - Aj^i{x)Aj^2{x)C{x) + Al^{x)E{x )], (4.11) 

^ Indeed, were any such curves to intersect at a;*, P would undergo a transcritical bifurcation along 
K“^(x*). Such bifurcations are not generically permitted by (Cl) and (C2). 
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where B, C and E are as in (I4.9|l . is called the j—fold bifurcation function of (j4.1|l . Let 
TZj be the zero-set of rj(a:) in We shall say that x* is a transverse intersection of 
and IZj if det A{x*) = rj(a:*) = 0 and 


rankH 


det y4(x*)A 

n(x -)) 


2 . 


Proposition 4.4 Let j G {1,2}. If x* is a transverse intersection of k{ 15 ) andIZj such 
that {Aj^i{x*) Aj^2{x*)) 7 ^ 0 and either 

(1) B{x*) = 0 and Aj,i{x*)C{x*) - Aj^2{x*)E{x*) = 0 or 

(2) B{x*) 7^ 0 and C{x*Y - 4 :B{x*)E{x*) > 0, 

then Pp undergoes a fold catastrophe at {x*, s*) for all s* such that 7p(s*) = ±ax*^p. 

Proof: By re-labeling the terms if necessary, we may assume Aj^i 7 ^ 0. There are then 
two possibilities. 

(i) If 5 = 0 and Aj^iC — Aj^ 2 E = 0, then 

= {X ; AiAaC + A^E = 0} = {A : A 2 = 0 or XiC + X 2 E = 0}. 


(a) If A 2 = 0, then L^* = {(Ai,0) : \iAj^i = 0} = {0} since Aj^i 7 ^ 0. But this 
contradicts the assumption dimLj,* = 1 . 

(b) If XiC + X 2 E = 0, then 


rank 




1 . 


(ii) If 5 7 ^ 0 and - ABE > 0, then 


In this case, 





-C ± VC2 - ABE 
25 



45r,- = (-2AI,-25 + Aj^iCf - - ABE) 

= (-2AI,-25 + Al ,-1 (C + ^C^-ABE^j 

■ (-2Al^-25 + Al ^-1 (C - VC^-ABE^^ = 0 

and so 

AI,-2 -C + VC^ - ABE A. 2 -C - VC'2 - ABE 

- — = - or- — = -. 

Alj-i 25 Alj-i 25 
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In either cases, L^* is contained in thus {0} 7 ^ La,* C K^.* and {0} 7 ^ L^* C K^*. 
As X* is a transverse intersection of fi;(3) and TZj, it is also a regular solution of (I4.8|l : 
(x*, La;*) then consists of fold bifurcation points of ilD, according to proposition 14.dl 
The desired conclusion follows from {±aa;*,p} = L^* = La,* n 7 p and from the correspon¬ 
dence between fixed points of Pp and ‘lines’ of hxed points of P. □ 

By construction, the bifurcation diagram of Pp is 27r—periodic in s. Consequently, 
elements of Cb must be (bounded) loops and elements of Coo must give rise to two 
cxD—catastrophes. Moreover, the number of fold catastrophes on any given C* G Cb can¬ 
not be odd as C could not be a loop were that the case. Finally, note that catastrophes 
cannot occur at 0 or as this would contradict (P2) and (P3). 

4-.3. The Bifurcation Diagrams 

Let Co, C^, So and s^ be as defined previously. Set rj G {0,.^}. By definition, 
goes through 77 at s = Sp. By (P3), the wedges’ angles he in (0, tt) or (0, tt] (when 
they exist), according to whether they record fold or 00 —catastrophes, respectively. Set 
ui = (ff -\- ipf and z /2 = 79 o -L Then Wq and overlap 

(i) in all four quadrants if and only if 1 ^ 1 , 1^2 > 7 r/ 2 ; 

(ii) in the first and third quadrants if and only if z/i > 7r/2 and 1^2 < 7 r/ 2 , and in the 
second and fourth quadrants if and only if < 7r/2 and z /2 > 7 r/ 2 . 

If Uj = 7 r/ 2 , the wedges do not overlap but their complement has zero measure in a 
neighbourhood of the origin. When the wedge angles do not exist, is a deleted 
neighbourhood of the origin in parameter space. 

4 . 3 . 1 . The case Co 7 ^ In this instance, it is sufficient to understand the bifurcation 
diagrams along a single curve: the full picture can then be obtained by combining 
the diagrams corresponding to Co and C^. When Cp G Cb, there are two (essentially) 
distinct generic possibilities. 

(i) If there is no fold catastrophe along C^, then the angles pfj do not exist and \Np is 
catastrophe-free deleted neighbourhood of the origin in parameter space. 

(ii) If there are 2k fold catastrophes along Cp, k > 0, then the angles are well-defined: 
Sri TTp are the s—values of the first fold catastrophes occurring respectively before 
and after rj along Cp. 

When Cp G Coo, there are two (essentially) distinct generic possibilities. 

(i) If there is no fold catastrophe along Cp, then the angles are well defined and 
Sp =F are the s—values of the cx)—catastrophes occurring respectively before and 
after rj via Cp. 
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Figure 3. Partial bifurcation diagrams of Pp when Cq ^ C^. Only one branch is 
shown. The square represents the origin or ^ and the circles and arrows indicate fold 
and oo—catastrophes, respectively. 


(ii) If there are k fold catastrophes along k > 0, then the angles are well- 
dehned: if all the fold catastrophes he on one side of p (say s > Srj) along then 
Sn =F are the s—values of the cxo—catastrophes occurring before p and the hrst 
fold catastrophe after p along respectively (or vice-versa). Otherwise, =F 
are the s—values of the hrst fold catastrophes occurring respectively before and 
after rj along Cr^. 

Some corresponding qualitative bifurcation diagrams are shown in hgure 3. 

4-3.2. The case Cq = In this instance, the bifurcation diagram must pass through 
0 at s = 0 and ^ at s = 11 / 2 . When Oq = G Cp, the number of fold catastrophes 
along the curve is even; there are then three (essentially) distinct generic possibilities. 

(i) If there is no fold catastrophe along Cq = then the angles and do not 
exist and Wq = are catastrophe-free deleted neighbourhoods of the origin in 
parameter space. 
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(ii) If there is an odd number of fold catastrophes between the origin and ^ along 
Co = C^, then the angles T^o and 7 r/ 2 =F 9 ?^ are well-dehned: they are the s—values 
of the hrst fold catastrophes occurring respectively before and after 0 and ^ via 
Co = Q. 

(hi) If there is an even number of fold catastrophes between 0 and ^ along Co = C^, the 
situation is much as described in ( 2 ), save for the fact that Cq = C^ is not a loop 
in the bifurcation diagram of Pp. 

When Co = G Coo, there are two (essentially) distinct generic possibilities. 

(i) If there is no fold catastrophe along Co = C^, then the angles (po and are well 

dehned and are the s—values of the cx)—catastrophes occurring respectively 

before and after 0 and ^ via Cq = C^. 

(ii) If there are k fold catastrophes along Cq = C^, A; > 0, then the angles ipo and 
are well-dehned: if no fold catastrophe lies between 0 and ^ along Co = C^ then 
=F<yi^o and 7 r /2 =f are determined as in the case C^ G Coo, item (2) (see d. If 
there are fold catastrophes between 0 and ^ along Co = C^, then ipQ and n/2 — 
are the s—values of the hrst fold catastrophes occurring respectively after 0 and 
before ^ along C^. 

Some corresponding qualitative bifurcation diagrams are shown in hgure 4. 


4.4- The General Mapping 

The mapping (ED is not the most general mapping satisfying (PI) —(P3); one should 
instead study maps of the form 

P{xi A) = x -|- 27r |^Ai.7^o(^, Ai) -I- AiA2iT(x, A) -I- \2G^{x, A 2 )], (4.12) 

where 7 ^ 0 G JPq, fT and are real analytic in their variables and the jacobians 
AAa;^o(0, Ai) and X 2 ) have the particular form prescribed by proposition 14.51 

which is analogous to proposition 14.11 


Proposition 4.5 //JF(0, Ai) = 0, ^^(^, A 2 ) = 0, and if 


DxiPo{0, Ai) 


a(Ai) -&(Ai)\ 
b{Xi) a(Ai) ) 


and DxQi{i,X 2 ) 


c(A 2) —d(A2)A 
d(A2) c(A2) 


where a,b,c,d : W ^ M. are continuous in their variables and a( 0 ),c( 0 ) 7 ^ 0 , then there 
exists ci;* > 0 such that the map defined by dm satisfies conditions (PI) —(P3). 

Dehne .4,: x —> M 2 (M) by 


A{x, A) 


-^o(^, Ai)-|-^J7’(a:, A) ^^/(x. A)-|-^^(x, A 2 ) ; 


(4.13) 
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Figure 4. Partial bifurcation diagrams of Pp when Cq = C'f ■ The squares represent 0 
and and the circles and arrows indicate fold and oo—catastrophes, respectively. The 
apparent self-intersection is an artifact of the projection on the ||x|| — s plane; it does 
not, in fact, occur. 


fixed points of are then in one-to-one correspondence with solutions of 

^(a:,A)-A = 0. (4.14) 

Set A) = det^(x, A). According to Taylor’s theorem, there are appropriate 
functions A'lo, A'oi such that '${x, A) = 0) -|- AiA'io(x, A) -|- X 2 Kqi{x, A). 

Let X be such that 0) = 0, det , 0) 7 ^ 0 and Al(f, 0) 7 ^ 0. Then, by 

the implicit function theorem, there is a neighbourhood QJ C of the origin and 
a unique analytic function X : such that X( 0 ) = x, i 5 ^(X(A),A) = 0 and 

rank^(A:(A), A) = 1 V A G 2J. 
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Define = {A G 2T : Al(X(A), A) ■ A = 0}. A simple rank argnment shows that is 
dehned via a single equation in two real variables, with a regular solution at the origin; 
consequently, as a manifold. Six is one-dimensional. 

Let Cx = kerAl(x, 0). Then, there is a small neighbourhood il C 5(0, cn*) of the 
origin in parameter space for which 

{(X(A), A) : A G il n is a deformation of {(x. A) : A G il fl ; 

both ‘curves’ can be parameterized by the same Xj, j = 1,2. The preceding discussion 
shows that the hxed points of are in one-to-one correspondence with the hxed 

points of the (already studied) truncated map 

Vt{x, X) = X + 2-11 [AiX'o(x, 0) A2^^(x, 0)]. (4.15) 

Fold bifurcations persist under small perturbations [35, 36]. Similarly, a generic 
unbounded curve remains unbounded under small perturbations. 

Indeed, in the real projective plane, an element of Coo meets the line at inhnity in 
two points. Generically, these two points are distinct and a small perturbation will not 
change that fact, i.e the perturbed curve is still an element of Coo- In the non-generic 
case where the two points at inhnity are equal, a small perturbation will either cause 
the points to separate or to vanish entirely (reminescent of a fold bifurcation of points 
at inhnity), i.e the perturbed curve either stays in Coo or becomes hnite. 

Thus, catastrophes generically persist: as a result, the bifurcation diagrams of 
(ITTl) and are (locally) topologically equivalent for small parameter values A. 

Consequently, (I4.12j) has wedge-like regions corresponding to the wedge regions 
of (jHH). 

Finally, note that since P is the ‘linearization’ of V at the origin with respect to A, 
the wedge regions of iH) provide tangential ‘cones’ for the corresponding wedge-like 
regions 211^ of ()4.12j) . 

5. Numerical Simulations and Examples 

In this section, we illustrate and interpret the results of the preceding sections through 
various examples. As such, the emphasis lies with qualitative observations rather 
than with precise numerical analysis. First, we study systems of PDF from a (naive) 
numerical perspective: we observe spiral anchoring, as well as hysteresis and homotopy 
of the spiral tip. Finally, we provide a few examples of mappings of the form (ITOD 
together with their zero-level sets and partial bifurcation diagrams. 

5.1. PDE, FESB and Semi-Flows 

In this section, we examine systems of partial differential equations giving rise to semi- 
flows satisfying the FESB equivariance described in section 2. 
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Figure 5. Anchoring in (EH) with perturbations as in 115.211 . The spiral tip paths 
are plotted in black, the anchored perturbed rotating wave is shown in gray, and the 
squares indicate the location of the perturbation centers. 


The computations are carried out on a two-dimensional square domain [—30, 30]^ with 
200 grid points to a side and time-step At = 0.005 and Neumann boundary condition, 
using a 5-point Laplacian and i) an explicit Runge-Kutta 2—stage method of order two 
in section 15.1.11 and ii) Matsui’s fourth-order Runge-Kutta code based on Barkley’s 
EZ-Spiral in section 15.1.21 Throughout, centers of anchoring are found via fast Fourier 
transforms of the tip data. 

5.1.1. Spiral anchoring Consider the following small perturbation of the FitzHugh- 
Nagumo equations: 

ut = -v) + (t)i +Au, 

(5.1) 

vt = g{u + (5 -'^v - 02), 

where 


4>j{x) = A/2cos(0.057r)0.12/(xi — Cij,X 2 — C 2 j), j = 1,2, (5.2) 

Ci,i = 9, C 2,1 = 0, Cl ,2 = -10, C 2,2 = 5a/3 and 

f{x) = exp (—0.00086 [xf + . 

Each gj{x), alone, breaks translational symmetry but preserves rotational symmetry 
about (9,0) (for j = 1) or (— 10,5v^) (for j = 2). Note that both perturbations are 
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uniformly bounded on and that they go to 0 as ||x|| —>■ cx). Under these conditions, 
the flow of (EH) near a normally hyperbolic rotating wave is equivalent to the flow of 
some center bundle equation dZH). Thus, if spiral waves anchor at all, they will generi- 
cally do so away from either perturbation center. This is conhrmed in hgure 5, in which 
the transients anchor at what would be an otherwise unremarkable location. 

We now present the results of simulations on a reaction-diffusion system with 4 TSB 
perturbations. Set 

ut=^-{u- -v)+(j)i + Au, 

vt = g{u + f3 - -fv + 02 ), 

where = 0.3, (3 = 0.6, 7 = 0.5, and where 0 i ,02 are inhomogeneous terms which 
depend on x G and are dehned by 

01 (x) = gi{x) + g 2 {x) = 0.12fi{xi - 9,X2) - 0 . 10 f 2 {xi + 1 ,X 2 - 10 ), 

h{x) = gsix) + gA{x) = -0.12/1 (xi 10,a;2 - 5 ^/ 3 ) -h 0.08/3(0:1 - 10,0:2 - 10), 

where = 0 . 12 , A 2 = - 0 . 10 , = - 0 . 12 , B 2 = 0.08, 

fj(x) = exp (aj(x^ + X 2 )) , j = 1, 2 , 3, 

oi = —0.00086, 02 = —0.0008 and 03 = —0.0009. 

Each gj{x), alone, breaks translational symmetry but preserves rotational symmetry 
about Cl = (9, 0) (for j = 1), C 2 = (—1,10) (for / = 2), C 3 = (—10, 5-\/3) (for j = 3) and 
C 4 = (10,10) (for j = 4). Note that the four perturbations are uniformly bounded on 
and that they go to 0 as ||a:|| —> cxo. As predicted, anchoring takes place away from 
the Cj, / = 1,... ,4. 

The transients in hgure 6 appear to hrst (hyperbolically) approach some manifold 
along which they travel to the anchored perturbed rotating wave; this will be the topic 
of an upcoming paper. 

5.1.2. Homotopy and hysteresis of rotating waves Following [ 8 ], dehne the modihed 
Oregonator 

Ut = ^(u-u^- ifv + 0)^) + Au, 

Vt = {u — v) + 0.6Au, 

where / = 1.4, g = 0.002, = 0.05 and 0 is an inhomogeneous term which depends on 

X e M^. When 0 = 0, (j5.4|l has full Euclidean symmetry. 
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Figure 6. Anchoring in the FitzHugh-Nagumo equations (15.31) . The spiral tip paths 
are plotted in black, the anchored perturbed rotating wave is shown in gray and the 
black squares indicate the location of the perturbation centers. 
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In the following simulations, 0 is the sum of two Gaussian bells; 


0(x) = O'! exp 


||(a:i,X2) - (15,15)||- 
Pi" 


+ 02 exp 


11 ( 3 : 1 , 3 : 2 ) - (18.75,15)|p 
P2" 


with «!, 02 , Pi, P 2 e M, and Pi, P 2 7 ^ 0. Each gj{x), alone, breaks translational 
symmetry but preserves rotational symmetry about Ci = (15,15) (for j = 1) or 
C 2 = (18.75,15) (for j = 2). Note that both perturbations are uniformly bounded 
on and that they go to 0 as || 3 :|| —> 00. 

When a = (ai,0) 7 ^ 0, ()5.4jl is §0(2)^—equivariant. Similarly, ()5.4jl is SO( 2 )c 2 — 
equivariant when a = ( 0 , 02 ) 7 ^ 0 and trivially equivariant when ai,a 2 7 ^ 0 . 


Set Pi = P 2 = I and p* = 0.01. Along the path a(r) = 7 i(r) = p* (cos(r), sin(r))"'" 
in parameter space, D undergoes a homotopy of perturbed rotating waves, whose tip 
paths deform continuously from a circle centered at ci, when r = 0 , to a circle centered 
at C 2 , when r = 7 r /2 (see hgure 7). 


Along the path a^r) = 72 (t) = ]^ 7 i('r), however, the homotopy is replaced by hystere¬ 
sis. As the parameters vary along the path, (Q has an anchored perturbed rotating 
wave whose tip path deforms continuously from a circle centered at Ci. At Mi G 72 , the 
rotating wave jumps (discontinuously) to another anchored perturbed rotating wave, 
whose tip path deforms continuously from a circle centered at C 2 . 
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Figure 7. Homotopy of the spiral tip path in 15.411 . The first spatial coordinates of the 
anchoring centers are plotted against r; compare this image with the first bifurcation 
diagram on p. [201 


Following 72 in the opposite direction leads to similar behaviour, this time with the 
discontinuous jump taking place at ^2 G 72 , as can be seen in hgure 8. Consequently, 
there must be a third unstable rotating wave (which escapes detection by direct means) 
appearing and then disappearing at Ki and ^ 2 , respectively, in saddle-node bifurcations. 

In a sense, both of these occurrences have been predicted by the analysis provided 
in section^ consult, for instance, the hrst and third bifurcation diagrams on p. 1201 

5.2. Wedges and Catastrophes 

In this hnal section, we provide a partial catalogue of (partial) bifurcation diagrams for 
mappings of the form 


P{x, X) = X + 271 [AiFo(a;) -b A 2 G$(a;)], 


(5.5) 


where A e ^ = (2,2)^, 



aij,bij,Cij,dij E M, i + j > 1, and /o and are continuous functions such that ()5.5j) 
satishes proposition 14.11 we can then use the visual criterion of section 14.21 to under- 
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Figure 8. Hysteresis of the spiral tip path in 15.411 . The hrst spatial coordinates of 
the anchoring centers are plotted against r. The question marks interpolate (roughly) 
the unstable rotating waves. Compare this image with the third bifurcation diagram 
on p. 1201 


stand the nature of the bifurcation diagram of the associated map Pq.oi- Furthermore, 
the wedge angles can be read directly from the bifurcation diagram. 

In the figures of this section, Cq and are shown in black or gray. Fold catastro¬ 
phes are indicated by circles, cx)—catastrophes by arrows and the squares mark both the 
origin and 

5.2.1. The Elowyn-Bonhomme map The Elowyn-Bonhomme (EB) map is obtained 
from (jSini) by setting fo = g^ = 1, opi = 62,0 = Co,2 = di,i = 0(2,0 = 1, 00,2 = -1 and 
all other coefficients to 0; the corresponding k(3) is shown in figure 8. In this instance, 
Co G Cb and Q G Coo- For the EB map, a;* = Let p = 0.01 < a;*. Using a pseudo¬ 
arc length continuation algorithm (see [37] for details), a partial bifurcation diagram of 
Pp (ignoring all fixed point branches but those through rj aX s = Sp, for rj G {0,^}) is 
built: the results can be seen in figure 9. 

There are 6 fold catastrophes: two along Co and four along C^. Their location can 
be recovered directly from and TZj, j = 1,2, with the help of proposition 14.41 in 
hgure 9, the six intersections that satisfy the appropriate hypotheses are marked with 
circles. Each corresponds to one of the six fold catastrophes observed in hgure 7. 

Furthermore, two cx3—catastrophes occur via C^. The interesting values for the EB 
map are compiled in table 1 above. We continue by providing examples that highlight 
the various possibilities. 
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Table 1. EB map catastrophes along Co and C^. 


Curve 

Type 

X* 

s* 

Wedge Angle 

Co 

Eold 

(1.2483,-0.1286)^ 

5.9809 

^0 

« 0.3023 

Co 

Eold 

(0.2269,-3.4760)^ 

0.2308 


w 0.2308 


Eold 

(0.3371,3.1473)'^ 

1.1020 

Pi 

« 0.4688 

Q 

Eold 

(2.2769,0.2982)^ 

2.1125 

Pt 

« 0.5417 

Q 

Fold 

(-3.2933,6.1024)^ 

1.1581 


n.a. 

Q 

Fold 

(5.6733,-1.2807)^ 

1.9267 


n.a. 

Q 

Infinity 

n.a. 

1.0172 


n.a. 

Q 

Infinity 

n.a. 

2.3562 


n.a. 



Figure 9. On the left: zero-level set k(3) for the EB map. On the right: fixed point 
branches of Pq.oi (the apparent self-intersection on Cq is due to a projection onto the 
||a;|| — s plane). Both Cq (the curve through 0) and (the curve through are shown 
in black. 


5.2.2. The first example 

02,0 = 0 Oi^i = 1 ao,2 = —1 (^2,0 = 1 (^1,1 = 0 6 o,2 = 0 

C2,0 = 0 Ci^i = 0 Co,2 = 1 d2p = 1 di^i = 1 do,2 = 0 

fo{x) = exp (-(x? -b xl)/lQ) g^{x) = exp (- ((xi - 2f + {x2 - 2)^) /14). 

See figure 11 for a portion of fi:(3) and a partial bifurcation diagram of Po.oi- By- 
construction, the zero-level set for this hrst example is exactly the zero-level set of the EB 
map; however, their bifurcation diagrams are not topologically equivalent (compare with 
hgure 9). In this instance, the wedge angles record fold catastrophes on Cq (black) G Cb 
and cx)—catastrophes on (gray) G Coo- Note further that this map provides an 
instance when the anchoring wedges overlap. 
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Figure 10. Intersections of the zero-level sets /t(3) (thick black lines) and TZi and TZ 2 
(thin gray lines) of the Elowyn-Bonhomme problem. The squares represent 0 and 
the points that satisfy the hypotheses of proposition 14.41 are marked with circles. The 
light gray region is part of the planar set for which C'(x)^ — 4:B{x)E{x) < 0. 




Figure 11. The first example: k(3) (left), bifurcation diagram of Pq.oi (right). 


5.2.3. The second example 


2 

02,0 - 5 

n - 43 

- 10 

00,2 — 

17 

2 

^ 16 

99 

Co,2 = 

48 

C2,0 - -y 

^1,1 - “To 

5 


b 

d 


_ 51 

2,0 10 


_ 53 

2,0 10 



91 

^0,2 — 

12 

'1,1 = 

10 

5 

f 

61 

do,2 = 

99 

1,1 ~ 

10 

10 


/o(x) = 1 g^{x) = 1 


See figure 12 for a portion of /t(3) and a partial bifurcation diagram of Pq.oi- la this 
instance, Cq (black), Q (gray) G Coo, and record fold catastrophes while 

records an cx)—catastrophe. 
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Figure 12. The second example: /t(3) (left), bifurcation diagram of Pq.oi (right). 



Figure 13. The third example: k(3) (left), bifurcation diagram of Pq.oi (right). 


5.2.4- The third example 


„ _ 53 

02,0 — “To 

„ — 22 
^2,0 - y 


Ol,l — 

<=u = f 


28 
' 5 


O0,2 
Co ,2 


= -9 


(^ 2,0 — 
d2,0 = 


_9 

5 

-9 


bri = 


41 

'lO 


g — 49 
01,1 - To 


ho,2 

do ,2 


= —^ 
10 

= -10 


fo{x) = 1 g^{x) = 1 


See figure 13 for a portion of fi:(3) and a partial bifurcation diagram of Po.oi- In this 
instance, Cq (black), (gray) € Cb, and the wedge angles all record fold catastrophes. 


5.2.5. The fourth example 


02,0 — 2 

n — 13 

“ 1,1 - T 

Oo ,2 — 

13 

10 

h — 16 

02,0 - -y 

7 ) — 21 

“ 1,1 - “To 

h — 3 
Oo ,2 - 5 

r - ^ 

C2,0 - - 5 

^ — 3 

“ 1,1 10 

Co ,2 = 

29 

10 

^2,0 = 

'^ 1,1 “ 10 

do ,2 = “1 



fo{x) 

= 1 

= 1 
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Figure 14. The fourth example: k(3) (left), bifurcation diagram of Pq.oi (right)- 




Figure 15. The fifth example: k(3) (left), bifurcation diagram of Po.oi (right). 


See figure 14 for a portion of k(3) and a partial bifurcation diagram of Pq.oi- la this 
instance, Cq = E Cb, the wedge angles all record fold catastrophes and the anchoring 
wedges overlap. 

5.2.6. The fifth example 


02,0 = 

59 

10 

01,1 — 

olco 

oo,2 — —7 

03,0 = 

38 

5 

02,1 — 

27 

5 

Ol,2 — 

59 

10 

Oo,3 = 

1 

5 

b2fi = 

37 

5 

&1,1 = 

51 

10 

h — 9 

00,2 - -5 

bsfl = 

63 

10 

b2,l = 

41 

10 

bl,2 = 

-4 

bo,3 = 

43 

10 

C2,0 = 

11 

5 

Cl,l = 

36 

5 

1 

o' 

0 

C3,0 = 

61 

10 

C2,l = 

2 

Cl,2 = 

4 

5 

Co,3 = 

8 

5 

<^2,0 = 

49 

10 

di,i = 

38 

5 

d - 27 

“0,2 - 10 

dsfl = 

79 

10 

d2,i = 

69 

10 

di,2 = 

29 

5 

do,3 = 

41 

10 





fo{x) 

= 1 


= 1 







See figure 15 for a portion of k(3) and a partial bifurcation diagram of Pq.oi. In this 
instance, Cq = E Coo and the wedge angles all record fold catastrophes. 
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Figure 16. Fixed point branches of T^o.oi (left); anchoring wedges Wg (gray) and Wq 
(black) for the Elowyn-Bonhomme maps (right). 


5.2.7. The Elowyn-Bonhomme map revisited Let 5,^ Fq and be as in the original 
EB map (see p. and consider the map V dehned via dm, with 


T ](a;, Ai) — Fq{x) + Ai 


— yXi + 9X2 


9 2 

5-^1 


—9xi — ^X2 — 9x1 + — 10x2 


41 rv. _ 57 ^2 

2 o ‘^ 1"''2 io ”''2 

^9^ 10^2 


^^(x, A 2 ) — G^{x) + A 2 


26xi + |X 2 + -jXf - 15 X 1 X 2 


f Xi - 6 X 2 + '^x'i 


28 ^2 


2™2 ' 
5-^2 

'T* _i_ 

2X1X2 “T K X2 


and 


J{x,X) = 


67 I 32 I 18 _ 33 2 _ ™ ™ _ 8 2 
10 + -5-2:1 + -5-X2 yX^ X1X2 5X2 


M _L 
10 5-^1 


46 „ I 5 ™2 I 79 „ ™ I 73 „2 
-X2 + 2^1 “T 102:1X2 + io3:2 ^ 


Then, according to proposition 14.51 V satishes (PI) —(P3). By the preceding discussion, 
the bifurcation diagrams of the EB map P is topologically equivalent to that of V when 
A is close enough to the origin. 

As an example, let 0.01 < a;* = and dehne Po.oi as the restriction of V to the 
circle 7 o.oi('S) = 0.01 (cos(s), sin(s)) in parameter space. Compare the diagram shown 
(on the left) in hgure 16 with the corresponding diagram of section IP).‘ 2 . 1 I 

Finally, the wedge regions Wq of the EB map and 2Ilo of the revisited EB map are 
shown (on the right) in hgure 16, illustrating the last remark of section 4. 
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